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Abstract. 

We update constraints on the Hubble function H{4>) during inflation, using the 
most recent cosmic microwave background (CMB) and large scale structure (LSS) 
data. Our main focus is on a comparison between various commonly used methods 
of calculating the primordial power spectrum via analytical approximations and the 
results obtained by integrating the exact equations numerically. In each case, we 
impose naive, minimally restrictive priors on the duration of inflation. We find that 
the choice of priors has an impact on the results: the bounds on inflationary parameters 
can vary by up to a factor two. Nevertheless, it should be noted that within the region 
allowed by the minimal prior of the exact method, the accuracy of the approximations is 
sufficient for current data. We caution however that a careless minimal implementation 
of the approximative methods allows models for which the assumptions behind the 
analytical approximations fail, and recommend using the exact numerical method for 
a self-consistent analysis of cosmological data. 
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1. Introduction 

Cosmic inflation is the simplest and most robust paradigm capable of providing self- 
consistent initial conditions to the Hot Big Bang scenario [H [21 El HI [6], as well as 
a mechanism for the quantum-gravitational generation of primordial scalar (density) 
perturbations and gravitational waves [3 El El EEl EH CGI [13]. The Fourier power 
spectrum V-ji(k) of the former is observed today in the cosmic microwave background 
(CMB) and the large scale structure (LSS). Vice versa, at present the CMB and the 
LSS provide the only quantifiable observables which can confirm or falsify inflationary 
predictions. That is why matching concrete inflationary models to observations has 
become one of the leading quests in cosmology. 

In the standard inflationary picture, the amplitude of perturbations for a given 
comoving Fourier mode k depends crucially on the dynamics of inflation around the 
time of Hubble exit for this mode. Each Hubble exit time is conveniently parameterised 
in terms of the number of e-folds N before inflation ends. The relation between k and 
N depends very much on the overall energy scale of inflation. The ensemble of modes 
observable in the CMB and in the quasi-linear part of the LSS power spectra corresponds 
to a range AN ~ 10 called "observable inflation" . The total duration of inflation is a 
priori unlimited, but the number of e-folds between the time at which the presently 
observable Universe became as large as the Hubble radius and the end of inflation can 
only vary in the approximate range 30 < N < 60, that will be called "relevant inflation" 
throughout this paper. 

The literature on inflation constraints is plethoric. For simplicity, a majority of 
papers are restricted to the case in which inflation is driven by a single scalar field (f> (an 
inflaton) with a canonical kinetic term, some potential V{4>) and minimal coupling to 
Einstein gravity (however, each of these assumptions can be relaxed and has already been 
studied separately). Traditional works are based on the definition of spectral parameters 
(amplitude, index, possibly running) for density perturbations and gravitational waves. 
In a first step, these parameters are fitted to the data; in a second step, one tries to infer 
the class of inflationary models compatible with derived bounds on spectral parameters. 

In the last years, many works have gone beyond this approach, recognising that 
the introduction of spectral parameters puts already a strong theoretical prior on 
the models, and is by no means a necessary step. It is more realistic and equally 
efficient to fit directly to the data the (more fundamental) parameters governing the 
dynamics of inflation and/or the inflaton potential. Within the main stream (standard 
single field inflation), recently published analyses fall in two categories which are 
both interesting and complementary: either one assumes a particular model based on 
a definite form for the inflaton potential throughout relevant inflation, and derives 
constraints on the free parameters of this potential (top-down approach); or one 
employs a generic parameterisation of the potential V(4>) or another function governing 
inflationary dynamics, e.g. H ((/)), and tries to reconstruct this function from the data 
(bottom-up approach). The second approach aims at avoiding theoretical priors as 
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much as possible, and concentrating on what the data exactly tells us, although no 
parameterisation can be completely general: sharp features are usually excluded ah initio 
(obviously, all possible features cannot be accounted for with a reasonable number of 
free parameters). Even within the bottom- up approach, a distinction can be established 
between conservative analyses reconstructing only the part of V(4>) corresponding to 
observable inflation; and more aggressive analyses in which the potential (or the function 
H ((/>)) is extrapolated till the end of inflation, and subject to a prior on the minimum 
duration of relevant inflation (e.g. iV > 30). This more aggressive method should 
imply varying many more parameters, since in this case the parameterisation should 
be accurate over 30 to 60 e-folds instead of just ~ 10. The fact of extrapolating is by 
itself an extra theoretical prior, since cosmological data tell us essentially nothing about 
the era between observable inflation and Nucleosynthesis: the end of inflation could 
be subject to multi-field dynamics, experience phase transitions, be split into several 
non-contiguous short inflationary stages, etc. 

Here we address only the most conservative approach, i.e., the reconstruction 
of the inflationary dynamics during observable inflation, with a minimal number of 
assumptions. After the publication of WMAP results, this approach was followed in 
Refs. [H [151 ESI EH HE]. It was stressed in [H [T^ □!] and 02] that the quantity 
primarily constrained by the data is H(<p): hence, this function is the one which 
should be parameterised in some way and fitted to the data. The knowledge of 
H(<p) uniquely defines the potential V(4>), and Ref. [18] presented the collection of 
potentials V((p) corresponding to the ensemble of functions H((f>) allowed by current 
WMAP and SDSS LRG data. The results are still plagued by degeneracies: since the 
energy scale of inflation is unknown, the data favors a parametric family of inflaton 
potentials rather than a precise shape within the observable window. However, a lot of 
improvement is expected from the next generation of CMB experiments, especially the 
planck satellite. The reconstruction of observable inflation will improve spectacularly 
if primordial gravitational waves are observed by planck or another experiment, in 
the form of polarised I?-modes. This would fix the tensor-over-scalar ratio r, and hence 
the energy scale of inflation. Even without I?-modes, the Planck data would provide 
r-dependent constraints on the inflation potential of unprecedented precision. In this 
perspective, it is worth comparing the details and merits of each reconstruction method. 

Choosing to concentrate on the reconstruction of H(4>) during observable inflation 
does not fix the method entirely, in particular as far as the computation of the primordial 



spectra is concerned. The authors of [T 
two different forms, while those of [18)1 



[T5l [T6] employed analytic approximations of 
wrote a module appended to CAMB [24| and 
CosmoMC [25], which derives the numerical spectra for each new set of inflationary 
parameters by numerically solving the exact equations. It is interesting to study whether 
the difference between these methods is relevant given the precision of current and future 

| Numerical spectrum computations were also employed in various complementary approaches to the 
problem of constraining inflation, based on definite potentials [19j |20l [2TJ [22] or on a frequentist analysis 
with extrapolation of the potential throughout relevant inflation |23j . 
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data. Beyond the issue of perturbations, different methods could also differ through 
different assumptions concerning the parameterisation of the background evolution, 
and the exact number of e-folds during which this parameterisation is (explicitly or 
implicitly) assumed to hold and be compatible with accelerated expansion. The goal 
of this paper is to compare in details these different techniques, and to see how each 
difference impacts the constraints obtained from current data. Although current inflaton 
potential reconstructions are still dominated by degeneracies, a careful understanding 
will be necessary before applying these methods to the highly precise data expected in 
the next years. 

In the next section we will briefly review the theory of inflationary perturbations 
and discuss the exact approach, as well as two commonly used approximative methods 
for the calculation of the primordial perturbation spectra. In section [3J we will present 
the results of an analysis of current data and demonstrate that the bounds using 
the approximative methods with a naive prior differ significantly from the constraints 
inferred with the exact method. We will track down the cause of these differences and 
compare the accuracy of the approximations in section H] before we conclude in section 

m 



2. Background and perturbations in single field inflation 

The observable spectra of density perturbations and gravitational waves are directly 
related to the evolution of the Hubble parameter H = a/ a as a function of in the 
neighbourhood of an arbitrary pivot value 0*. The function H((p — 0*) can in principle 
be reconstructed from the data without any need to assume an explicit value of 0*. 
Each H((j) — 0*) defines a unique set {V(0 — 0*), ini } through 

S2-7T 2 1 2"7T 

- - <f>*) = TO - - - <M> 

M 2 

' ] H'io o.). (2.2) 



47T 

whenever ^ [§| (the prime denotes a derivative with respect to 0, and we have set 
GMp Y = h = c = 1). In Ref. [H], the defining quantity if (0 — 0*) was Taylor-expanded 
up to the cubic term: 

F(0 - 0,) =H m + #:(0 - 0,) + ^;'(0 - 0*) 2 + J#f(0 - 0*) 3 , (2.3) 

2 o 

which is equivalent to keeping the first three slow-roll parameters 

(2.4) 



47T 



H n ' 2 



H 



§ Such a singularity is never reached as long as H{(f)) is used as the defining quantity and has an 
analytic expression over the range considered. As mentioned in Ref. [18] , (f> = can be reached for a 
field value <pi on ly if H has a non-analytical expression like (H — Hi) cx (<fi — <f>i) 3 ^ 2 in the vicinity of 
(f>i. This cannot happen with the parametrisations used in this work (polynomial expressions for either 
H(4>) or H 2 ((f>)). In addition, reaching = would imply that the field changes direction with an 
opposite sign for H' , which is not compatible with the assumption of a single- valued function H(4>). 
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Ml H" , . 

'=1F¥' (2 ' 5) 

* = 16^ (2 ' 6) 
in the Hubble flow hierarchy [261 127] . as in the slow-roll reconstruction approach of 
Refs. [281 [HI [151 US]. Note that the universe expansion remains accelerated as long 
as e < 1. For practical purposes, any parameterisation of H(4>) could be used when 
fitting the data with a Bayesian MCMC analysis. Besides the issue of priors on 
inflationary parameters, each parameterisation corresponds to a different ensemble of 
possible inflationary models. One can wonder how much the final results (i.e., the range 
of allowed potentials) depends on the parameterisation. In the following analysis we will 
compare the results obtained by Taylor-expanding either H(<f)) or H 2 ((p) at the same 
order, using in both cases the same flat priors on the first slow-roll parameters expressed 
at the pivot scale 0*. 

Once the ensemble of possible inflationary models has been specified, the analysis 
still depends on the way to calculate the perturbation spectra in single-field inflation, 
and on a theoretical prior on the duration of inflation. In this section we will give a 
very brief summary of three different approaches used in Refs. [HI [151 ESI EE] . Keep in 
mind that throughout this paper, we are working under the assumption of a minimal 
prior. In other words, we do not impose any lower bounds on the number of e-folds 
of inflation, we only demand that the spectra of a model can be calculated with the 
respective methods. So, besides possible differences in the accuracy of the resulting 
spectra, the default implementation of the three methods will also differ in the range of 
parameter values that would be excluded straight away. 

2.1. Exact spectra via mode equation 

Given H as a function of (ft during inflation, the spectrum of curvature perturbations Vn 
and gravitational waves Vt can be calculated exactly by integrating the scalar/tensor 
mode equation (see, e.g,. |29j): 

d 2 £s,T 



drf 



+ 



k 2 1 d 2 z s ,T 

z s ,t drf 



£s,t = (2.7) 



with rj = J dt/a(t) and Zs = a<p/H for scalars, Zt = a for tensors. The evolution of the 
background is determined by 

Ml dH , . 

47T d(j> 

The £s,t are usually taken to be in the Bunch-Davies vacuum when they are well within 
the horizon, and their evolution needs to be tracked until |£s,t|/2s,t converges to a 
constant value, in order to define the observable spectra: 

*** ( , 9) 



2tt 2 z% /c ' ttM^j 
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In principle, observable inflation could be interrupted for a very short amount of time, 
resulting in characteristic features in the spectra. In the mainstream approach, this 
situation is not considered for simplicity. Actually, the numerical module used in 
Ref. [18] eliminates models violating e < 1 at any point during the period of time 
the mode equation is integrated. More precisely, for any wavenumber in the range 
[&mm, &max] = [3 x 10 _6 , 1.2] Mpc -1 needed by CAMB (the pivot scale being fixed at 
fc* = 0.01 Mpc -1 ), the module integrates Eq. (12 .7p from the time at which k/aH = 50 
and until [d In "P^r/d In a] < 10~ 3 . If for a given function H((p — 0*) the product aH 
does not grow monotonically by a sufficient amount for fulfilling the above conditions, 
the model is rejected (we recall that it is equivalent to impose that aH grows or that e 
is greater than one). 

The condition [d In P^/d In a] < 10~ 3 is motivated by our desire to obtain 
a 0.1% accuracy in the power spectra. The error made on Vk,t by stopping the 
integration of perturbations at a finite time can be estimated analytically, comparing the 
amplitude of the decaying mode to that of the non-decaying mode for TZ or gravitational 
waves during inflation. The decaying over non-decaying mode ratio evolves in a first 
approximation like a -1 , i.e., like e~ N . Hence, a few lines of algebra show that the 
derivative [dlnP/dlna] is a good approximation for the relative error [AV/V] produced 
by stopping integration at a finite time. Other parameters governing the precision of the 
power spectra calculation (like the step of integration, the choice of the initial integration 
time for each mode, etc.) where chosen in such way that the above source of error is 
the dominant one. 

The numerical evaluation of the spectrum involves solving equations (12.71) for each 
value of k, but this does not increase the total running time of a Boltzmann code like 
CAMB by a noticeable amount. Nevertheless, there exist a number of approximations in 
the literature, which simplify the calculation considerably. 

2.2. Approximation I 

This method was employed in [THUS], and relies on the validity of the analytical slow-roll 
approximations, 



with C\ = —2 + In 2 + 7, where 7 is the Euler-Mascheroni constant. These equations 
were first derived in [30] and are accurate only to first order in the slow-roll parameters, 
assuming additionally that e and n are constant. Here, one only needs to solve one 
differential equation to determine <f)(k), 




(2.10) 



(2.11) 



din A; 



d<f) 



Mpi 

2^ 1 - e 



(2.12) 
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assuming = 0*. Once <fi{k) is known, the slow-roll parameters and hence the 

spectrum can be evaluated for each value k. In this approach, the evolution of Eq. (12.121) 
has to be followed throughout the observable range of wavelengths. If e > 1, equation 
(I2.12p will diverge, so models with e > 1 within this range will have to be excluded 
when using this method. If however, the inflationary condition were violated just before 
or after this range, the model would not be ruled out, and the resulting spectra would 
likely be inaccurate. 

2.3. Approximation II 

This method is based on the usual Taylor-expansion of the spectra in log-space around 
a pivot scale fc* (see e.g., [3T1 132] ). 

\nV n ~ In A s + (n s - 1) In (k/K) + ^a s (In (k/K)) 2 , (2.13) 

\nV T — In A^ + nj< In (k/k*) , (2.14) 

with the spectral indexes ng/x, and the running of the scalar tilt as given by their 
second-order slow-roll expressions 

n s ~ 1 + 277 - 4e - 2(1 + C 2 )e 2 - i(3 - 5C 2 )e V + 1(3 - C 2 )£, (2.15) 

as * -j^ + 8e 2 - 10e V + + ^ ^) , (2.16) 

n T ^ -2e - (3 + C 2 )t 2 + (1 + C 2 )er], (2.17) 

where C 2 = 4(ln2 + 7) — 5. The slow-roll parameters only need to be evaluated at a 
field value 0*, corresponding to the time when k* leaves the horizon. A$ and A^ are 
calculated from equations (12.101) and f)2.1ip . and the spectra follow directly. One does 
not need to solve any differential equations here, so the numerical implementation of this 
method is by far the simplest of the three. However, due to the additional assumption 
on the shape of the spectrum, it becomes increasingly inaccurate the further one goes 
away from the pivot scale. 

Apart from that, in the spirit of choosing a minimal prior one would typically 
rule out only those models that break the e < 1 condition at the pivot scale, thus 
allowing regions in parameter space in which inflation would break down even within 
the observable range and making the prediction of the spectra for these models extremely 
unreliable. 

3. Constraints from current data 

In this section we present the constraints on inflationary parameter space from a selection 
of current observations, comprising CMB data from the WMAP [33J, [M] , Boomerang 
[35| l36| 137] and ACBAR [38] experiments, complemented by the galaxy power spectrum 
constructed from the luminous red galaxy sample of the Sloan Digital Sky Survey [39|. 
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Table 1. Minimal 95%-credible intervals for the slow-roll parameters in the H- and 
if 2 -expansion schemes, using the exact method for calculating the spectra. 







H 




ff 2 


e 





0.028 





-> 0.023 


V 


-0.035 


-> 0.046 


-0.035 


-► 0.039 




-0.0026 


-> 0.028 


-0.0053 


-► 0.027 



We analytically marginalise over the luminous to dark matter bias b 2 and the nonlinear 
correction parameter Q n \. 

We consider a ACDM-model with eight free parameters, on which we impose flat 
priors. Four of these parameters determine the initial perturbation spectra: the scalar 
normalisation In [10 10 v4s], and the first three slow-roll parameters: e, rj and £, evaluated 
at the pivot scale fc* = 0.01 Mpc -1 . We emphasise once more that the numerical 
computation of perturbations does not refer to any slow-roll expansion, and remains self- 
consistent even when the field is not rolling very slowly. The fact of varying parameters 
which coincide with the usual slow-roll parameters is just a choice of prior in parameter 
space, which is particularly convenient for two reasons: first, the posterior is well- 
behaved with respect to these parameters and the convergence of the chains is achieved 
in a reasonable amount of time; second, it facilitates comparison with other works. 
The remaining four parameters are the baryon density u>b, the cold dark matter density 
co>dm, the ratio of sound horizon to angular diameter distance at decoupling 9 S , and 
the optical depth to reionisation r. We use a modified version of the Markov-Chain- 
Monte-Carlo code CosmoMC [2H [25] to infer constraints on the free parameters of the 
model. The inflation module was made publicly available by the authors of Ref. [H] at 
http : //wwwlapp . in2p3 . f r/~valkenbu/ inf lationH/ 



3.1. Expansion in H vs. expansion in H 2 

We first check the impact of changing the parameterisation of H(<p) (i.e., the precise 
ensemble of inflationary models considered) from a Taylor-expansion of order 3 in 
H(<p) to the same expansion in H 2 ((p). In both cases, we used the same priors on 
inflationary parameters: hence the difference only resides in the fact that slightly 
different background evolutions can be achieved in both cases. The differences are 
summarised in table Q] and turn out to be very minor. This preliminary analysis shows 
that the parametric form assumed for H((ft) within the observable window has a minor 
impact. Significant differences could only be expected if the choice of parameterisation 
of H(4>) would allow much more freedom in one case than in the other. 

In the remaining part of the paper, we shall therefore stick to the Taylor-expansion 
in H((j)) and perform three independent analyses, calculating the primordial spectrum 
either by exactly solving the mode equations, or using one of the two approximations 



How to constrain inflationary parameter space with minimal priors 



9 




n £ 



Figure 1. This plot shows the one-dimensional marginalised posterior distributions 
for the free parameters of the model. The black lines represent the results of the exact 
solution of the mode equation, red lines are approximation I and purple (dashed) lines 
correspond to approximation 77. 

discussed in sections 12.21 and 12. 31 
3.2. Approximations vs. exact spectra 

Our results are presented in figures CD and [3 We do not find any significant differences 
in the posterior probabilities of r, 9 S , uj h and tUdm- The four parameters that determine 
the primordial spectra, however, are more sensitive to the method used. Note that the 
exact method produces tighter bounds on the slow-roll parameters, particularly on £. 

This also has important consequences on the inferred values of derived 
phenomenological parameters, such as the spectral index and its running. As can be 
seen from table [2J the exact method yields significantly stronger constraints on these 
two parameters. 

4. Why the difference? 

There are potentially two reasons for these observed discrepancies. The first one is 
that the accuracy of the approximations could be insufficient within their respective 
"allowed" parameter space and lead to a serious bias in the parameter estimates. Note 




Figure 2. This plot shows the 68%- and 95% credible regions of the two-dimensional 
marginalised posterior in the (e, 77)- (top left), (e, £)- (t°P right), and (77, £)-planes 
(bottom). The black lines denote the results of the exact solution of the mode equation, 
red lines are approximation / and thin lines correspond to approximation II. 



Table 2. Minimal 95%-credible intervals for the spectral index and the running at a 
scale of k* = 0.01 Mpc -1 . Note that these are derived parameters and the results are 
not independent of the choice of pivot scale. 



exact approximation I approximation II 

n s 0.959 -> 1.049 0.960 -> 1.078 0.960 -► 1.087 
a s -0.063^ 0.001 -0.084^0.009 -0.098^0.003 



that the discrepancy occurs mostly in regions of parameter space where £ is large. The 
larger £, the more one would expect the accuracy of the approximations to degrade. 
However, given that the approximations are expected to be accurate to order £, i.e., not 
worse than ~ 10%, an effect as large as the one we observe seems rather unlikely. 

The second reason is slightly more subtle: as we discussed in section [2j the three 
methods differ in their implicit prior on the space of models. While approximation II 
requires e < 1 only at the pivot scale, approximation I needs us to demand that this 
condition be fulfilled in the entire observable window of ~ 10 e-folds, corresponding to 
Hubble exit for modes in the [k min , fc max ] range. In the exact numerical approach, we 



How to constrain inflationary parameter space with minimal priors 



11 



require e < 1 for the whole integration time, which starts when aH/k m i n =1/50 instead 
of one, and ends when |^s,t|/^s,t freezes out, i.e., a few e-foldings after aH/k = 1, 
corresponding to an even more restrictive prior. 

It was pointed out in Refs. [151 HO] that for models with large positive values of 
£ (> 0.05) and no higher derivatives, inflation tends to end within a few e-foldings of 
the pivot scale leaving the horizorjjj]. This is consistent with our results, since the more 
restrictive priors lead to tighter bounds on £. 

4-1. The prior issue 

To verify that the differences actually stem from the choice of priors and not from a 
lack of accuracy, we post-processed our Markov chains of the approximate methods, 
discarding all models for which inflation is interrupted in the range of wavelengths 
required for the exact calculation. 

In a first step, we remove only those models for which the inflationary condition is 
violated before the pivot scale leaves the horizon, when aH is in the range [fc m j n /50, A;*]. 
These are models for which the assumption of the Bunch-Davies vacuum initial condition 
is violated at least for the largest observable wavelengths. Only a mere 0.02% of the 
models in the chains using approximation I, and 0.01% for approximation II, fall victim 
to the cuft This is probably connected to the dislike of the data for models with large 
negative £, which is required if we want inflation to start only just before the observable 
range. 

Imposing the same additional prior as in the exact method (that inflation holds 
till the time of freeze out for each mode), ~ 20% of the approximation / points and 
~ 34% of the approximation // points are removed. After weeding out the bad models, 
the bounds of the approximations perfectly agree with the ones derived using the exact 
method, their marginalised posteriors are virtually indistinguishable. This confirms 
our suspicion that the different priors are responsible for the discrepancy between the 
methods. 

4-2. Comparison of accuracy 

Having seen that the prior plays a very important role, it is nonetheless interesting to 
take a closer look at how the approximations compare to the exact method in terms of 
accuracy. 

In order to compare the accuracy of the three methods described previously, we took 
the 95% best-fitting spectra obtained using approximations I and II and for each model 
in the chains we again computed the curvature spectrum in either of the approximations 
and numerically in order to compare. We searched for the maximum discrepancy 

|| If higher derivatives are present this conclusion can be weakened, see, e.g., pHl W2\ . 
% Here, and in the following, we quote a weighted fraction of models, i.e., ■ w° &d ) / ^ Wi), where Wi 
are the statistical weights of the points in the Markov chains, and vj° ad are the weights of the models 
killed by the prior. 
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Figure 3. Left: Curvature spectrum obtained from the exact numerical method 
or from approximation 7 or 77, for the most discrepant models in the range 
k 6 [3 x 10~ 6 ,1.2] Mpc -1 selected among the 95% best-fitting spectra. Right: The 
corresponding relative differences. The discrepancy is large only for scales close to 
&min, which have a relatively small statistical weight. 



Table 3. Value of the slow-roll parameters at the pivot scale for the models referred 
to as A and B in the text. The corresponding curvature spectra computed with each 
method are shown in Figure [3] 



model e rj £ 

~k 1.1446 x 10" 3 -4.1728 x 10" 3 3.4911 x 10" 2 
B 2.6159 x 10" 2 7.5254 x 1(T 2 3.3839 x 1(T 2 



between the approximated spectrum (with method lor IT) and the numerical one, in each 
of the two ranges [k min , h*\ and [k # , k max \, with k min = 3 x 1CT 6 , k* = 0.01, fc max = 1.2, in 
units of 1/Mpc, corresponding to the range of wavelengths the data are most sensitive 
to. Note that the spectrum for a model can only be computed numerically if the model 
meets the prior condition on the duration of inflation. Hence the comparison done here 
is for models that are already preselected by that particular prior, whereas in the actual 
chains many points exist that give a much larger discrepancy due to the different prior. 

In Figure [3] we plot the most discrepant models in the [k*, k m3iX \ range, that we 
call A (for approximation T) and B (for approximation IT). The corresponding slow-roll 
parameters evaluated at the pivot scale are given in Table [3l For approximation 77 we 
find a maximum discrepancy in Vn of 83% below k*, and 19% above (model B). For 7~, 
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the difference reduces to 33% below fc* and 8% above (model A). So, approximation / 
is doing better on both sides of the pivot scale. 

Indeed, it appears that approximation / captures the spectrum shape very well, but 
not its amplitude. This is not surprising since in i, the expression for the amplitude is 
first-order in the slow-roll expansion. The next-order contribution to this approximation 
should include the parameter £, with a positive coefficient. Since the data allows large 
values of |£| only when £ is positive, approximation /yields a systematic underestimation 
of the amplitude. This approach could be straightforwardly improved by computing the 
curvature amplitude at the next order. In contrast, approximation //tends to give the 
wrong shape, but since it includes one more order in the slow-roll expansion it can give 
a better estimate of the amplitude for k not too far from k*. 

However, the main result of this section is that the difference between the various 
spectra is very small, since a large discrepancy is only encountered on scales close to k min , 
which have a relatively small statistical weight in the process of accepting or rejecting 
a model. For both approximations the error in the larger part of the spectrum is of the 
order of 5%. Current data does not reach such good sensitivity, especially if we keep 
in mind that for the largest k values the curvature spectrum is mainly constrained by 
the SDSS data, which is always marginalised over an unknown bias parameter. For the 
particular models shown here, the difference in the effective x 2 obtained when fitting 
either the approximated or the numerical spectrum to the data is |A%g ff | = 6.7 for 
approximation / (model A) and |AXg ff | = 5.9 for approximation // (model B). One 
should keep in mind though, that these are just the most extreme deviations, at the 
edge of allowed parameter space. On average, the inaccuracies are too small to have a 
significant effect on the inferred bounds. However, with future datasets one can expect 
the |Axg ff | to become even larger, possibly resulting in biased estimates. Hence, we 
recommend using the exact numerical approach, since it does not make the analysis 
longer or more difficult. 

The conclusions reached in this section apply to a particular class of inflationary 
models, namely those described by Eq. (12.31) with parameter ranges limited by current 
WMAP and SDSS results. Allowing for more freedom in H(4>), one would expect 
stronger deviations between the analytical and numerical approaches. Conversely, 
imposing a constraint on the total number of e-folds for relevant inflation, one would 
select models which are deeper within slow-roll and obtain even smaller discrepancies. 
We have limited our discussion to the scalar spectrum, since there is, at present, no 
evidence for anything but a subdominant tensor contribution in the available data. 
This may of course change if a primordial 5-mode polarisation of the CMB is detected 
in the future. 

5. Discussion 

We have compared various alternative methods for putting constraints on the observable 
window of inflationary dynamics, assuming single-field inflation with a smooth 
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behaviour. One could fear that the results would depend very much on the way to 
compute the spectrum, or on the parameterisation of H(<p) (i.e., the ensemble of models 
considered). We point out that with current data these differences are subdominant. 
The results are mainly affected by the exact prior on the background. By focusing on 
the allowed window, one hopes to be conservative and to get results dictated by the data 
only; however these results are very sensitive to the edges of the interval in field space 
over which accelerated expansion is required. In other words, the upper bound on £ is 
given by this constraint rather than the shape of the perturbation spectrum. Current 
analyses agree with each other, but only within a factor two, due to this difference. This 
dependence of the final results on the choice of prior may sound worrisome, particularly 
if it is one's aim to keep the analysis as general as possible. 

It seems as though the priors we initially chose to use with approximations / and II 
are less restrictive, and would therefore lead to a more conservative result. Unfortunately 
however, combining these methods with their respective priors on the space of allowed 
models is fraught with a severe consistency problem. Both priors allow models in regions 
of parameter space where the approximations are known to break down, and yield results 
for the spectrum that cannot be trusted. Approximation II, for instance, allows models 
in which inflation is interrupted within the observable range. Such an event would lead 
to very distinct signatures, like a cutoff, yet the approximation would still predict a 
smooth spectrum. It is therefore sensible to expect that inflation lasted at least over 
the whole observable range. But even that will not be sufficient: if inflation started 
only just before the observable range, the assumption that the modes start out in the 
Bunch-Davies vacuum can no longer be justified, and the approximations fail. If, on 
the other hand, inflation ends just after the smallest observable scale leaves the horizon, 
the corresponding mode will not have time to freeze out. In fact it would re-enter the 
horizon right away, and there is no telling (without making further assumptions) in what 
shape the spectrum would arrive at later times when it is relevant for the determination 
of the CMB anisotropy spectra. 

Hence, it is reasonable to demand a proper vacuum initial condition and a freeze- 
out of the modes. One could in principle further limit the space of allowed models by 
constraining the minimum number of e-folds before the end of inflation to a certain 
number, usually taken to be > 30. However, this would require a daring extrapolation 
of our simple Taylor-expansion over a huge range of e-folds, where even a tiny higher 
derivative of the Hubble parameter would eventually take over. In this paper, therefore, 
we did not want to make any additional assumptions about what happens after the 
freeze-out. 

There are, however, two points at which the prior of the self-consistent numerical 
approach is slightly arbitrary, corresponding to the two end points of the interval over 
which we track the background dynamics. The first one is the large scale end, determined 
by the time at which we choose the initial conditions for the k m [ n . We picked k min /50 as 
a starting point, but other choices may have been equally good. Fortunately, the data 
conspire to make this choice have little impact on the final results: only about 0.1% of 
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the models in our chains generated with the laxer prior can be rejected due to inflation 
starting "too late". 

The more critical issue is the small scale end of the interval. Its choice is connected 
with the question when a mode can be considered to have frozen out, which is set by the 
limiting value of [dlnP^y/dlna] at which we stop integrating. The final results for the 
posterior are mildly dependent on the choice of this limit, which should eventually be 
chosen such that the resultant uncertainty in the spectra is smaller than the sensitivity 
of the data. For this reason, we choose here [d In P^r/d In a] < 1CT 3 , corresponding to 
a 0.1% accuracy in the power spectra. 

We would like to emphasize once again that the differences in the results are not 
inherent to the approximations used, but rather due to the attempt to implement them 
with a minimal prior. Our results also show that if one were to impose a non-minimal 
prior on the number of e-foldings beyond the observable range, the approximations 
would lead to the same results as the exact method. 

In the future, we expect more robust constraints from high-precision experiments, 
such as, e.g., the Planck satellite. In turn the difference between the various methods 
for computing the spectra will become even more relevant. In the light of our results, 
we recommend using the exact numerical approach for a self-consistent analysis of 
inflationary dynamics. 
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